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ABSTRACT 

We present a model for the formation of relativistic vortices (strings) at a quench, 
and calculate their density and correlations. The significance of these to early 
universe and condensed-matter physics is discussed. 



1. Introduction 

Despite our improved understanding of the very early universe the nature of 
the phase transitions in the Grand Unification era is largely unknown. This would 
seem to be a serious drawback in our modelling of the universe since, in all but the 
baldest inflationary scenarios, an imprint of these transitions is visible today. In 
particular, it has been argued by Kibbletl and others@ that the large-scale structure 
of the universe can be attributed to cosmic strings (vortices in the fields) formed at 
that time. 

However, it has been suggested that the initial conditions of any string network 
are largely washed out after a few expansion times, at which the network is assumed 
to approach a scaling regime with a few large loops and long strings per horizon 
volume continuing to produce smaller loops by self and mutual intersection. This 
has seemed to obviate the need for any detailed description of the early microscopic 
dynamics that set the boundary conditions for the latter classical picture. 

Despite this, there are two circumstances in which this lack of detailed initial 
information leaves us at a loss. The first concerns the density of the defects formed. 
Although not of immediate importance to the cosmologists, for the reasons given 
above, it does relate directly to the interesting proposition! that the production of 
vortices in superfluid A He and 3 He may share many attributes with the produc- 
tion of cosmic strings. To date, the primary data from such experiments concern 
defect densities!!!. If there is a similarity between the production mechanisms in 
these different regimes of high T (low chemical potential) and low T (high chemical 
potential) we should be able to predict the density of vortices in both superfluid 
experiments and the early universe. 
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Secondly, and of direct importance to the astrophysicists, the scaling solutions 
for the early universe mentioned above require the presence of some 'infinite ' string, 

1. e., vortices that do not self-intersect. The presence of such string is, in part, 
determined by the initial conditions and it is important to know whether it is 
present for reasonable models. For example, it has been suggested that vortices 
produced by bubble nucleation in a strong first-order transition will only form small 
loopsH The tendency for string to form loops should be visible in string density 
correlation functions. In a condensed matter context, the same correlation functions 
should enable us to estimate the superflow that would occur at a superfluid quench 
from fluctuations alonea. 

In this talk I shall show how these problems can be addressed in a model of 
vortex formation by unstable long-wavelength Gaussian fluctuations at a tempera- 
ture 'quench'. For simplicity, flat spacetime is assumed. Relying on a slow 'rollover', 
this excursionary model can only describe weak coupling systems for the short times 
while the domains are growing before the defects freeze out. Although this is unsat- 
isfactory for most early universe applications and for low-temperature many-body 
systems, we know in principled how to include back-reaction (still within the con- 
text of a Gaussian approximation) to slow down domain growth prior to the field 
fluctuations spreading to the ground-state manifold. If the results of Ref.7 are a 
reliable guide to our model then the conclusions that we shall draw are likely to 
survive. Greater detail is given in recent .work by myself and my collaborators, Tim 
Evans, Alasdair Gill and Glykeria Karrap'B, from which this talk is drawn. 

2. Vortex Distributions 

Consider the theory of a complex scalar field 0(x, t). The complex order param- 
eter of the theory is (0) = f]e ta and the theory possesses a global 0(2) symmetry 
that we take to be broken at its (continuous) phase transition. 

Initially, we take the system to be in the symmetry-unbroken (disordered) phase, 
in which the field is distributed about = with zero mean. We assume that, at 
some time t = to, the 0(2) symmetry of the ground-state (vacuum) is broken by a 
rapid change in the environment inducing an explicit time-dependence in the field 
parameters. Once this quench is completed the 0-field potential V((f>) = — M 2 |0| 2 + 
A|0| 4 is taken to have the familiar symmetry-broken form with M 2 > 0. 

In practice we expect that, as the complex scalar field begins to fall from the 
false ground-state into the true ground-state, different points on the ground-state 
manifold (the circle S 1 , labelled by the phase a of (0)) will be chosen at each point 
in space. If this is so then continuity and single valuedness will sometimes force the 
field to remain in the false ground-state at = . For example, the phase of the 
field may change by an integer multiple of 2n on going round a loop in space. This 
requires at least one zero of the field within the loop, each of which has topological 
stability and characterises a vortex (or string). As to the density of the strings, 
if the phase a is correlated over a distance £, then the density of strings passing 



through any surface will be 0(£~ 2 ) i.e. a fraction of a string per unit correlation 
area. On the completion of the transition a network of strings survives whose further 
evolution is determined by classical considerations as the field gradients adjust to 
minimise the energy. 

The question is, how can we infer these late-time string densities and the density 
correlations from the microscopic field dynamics? The answer lies in the fact, noted 
earlier, that the string core is a line of zeroes of the fields $ a (a=l,2). This is 
equally true for both relativistic and non-relativistic 0(2) theories. The problem is 
solved if we can identify those zeroes which will freeze out to define the late-time 
vortices. This will require careful winnowing, since it is apparent that quantum 
fluctuations lead to zeroes of the fields on all distance scales (even in the disordered 
phase). However, most of these zeroes will be transient. Two levels of screening 
are required before the relevant zeroes can be identified. Firstly, those zeroes whose 
positions vary rapidly, or which annihilate one another on time scales of 0{mT v ), 
where m provides the mass-scale, can be discounted. Secondly, since m~ l also 
measures the typical vortex width, zeroes on scales smaller than this should also be 
ignored. For the moment we ignore this problem, and count every zero. 

To see how to proceed, consider an ensemble of systems evolving from one of a 
set of disordered states whose relative probabilities are known, to an ordered state 
as indicated above. As a prologue to the problem in hand, counting the zero line 
densities appropriate to vortices, we consider the much simpler case of counting the 
zeroes of a real field <&(#) in one t space dimension. On implementing the change 
of state on the line, these zeroes fluctuate and annihilate, but some of them will 
come to define the position of 'kinks' (at which <&'(x) > 0), some the position of 
'antikinks' (at which < 0), the one-dimensional counterparts of vortices and 

'anti'-vortices. 

Suppose, at a given time, the zeroes of 3>(x) occur at x — xi,X2, ■■■■ It is useful 
to define two densities. The first, 



is the total density of zeroes, not distinguishing between kinks and antikinks (by 
which we now mean zeroes at which the field has positive or negative derivative). 
The second is the topological density, 



where = sign($' (xi)) , measuring (net) topological charge, the number of kinks 
minus the number of antikinks. 

^Because of the peculiarities of one spatial dimension, that would confuse the issue, we pretend that 
we are examining a one-dimensional subset of a real field in higher dimensions and ignoring other 
degrees of freedom. 





(2) 



Equivalently, in terms of the $-field, the totai density is 



p(x) = 6Mx)}\&(x)l (3) 

since <&'(x) is the Jacobian of the transformation from zeroes to fields. Similarly, 
the topological density is 

p{x) = 6[$(x)]&( x )- (4) 



Analytically, it is not possible to keep track of individual transitions, but we 
can construct ensemble averages. If the phase change begins at time to then, for 
t > t , it is possible in principle to calculate the probability Pt[<&) that <fi(x,t) (the 
counterpart of <f>(x,t)) takes the value at time t. Ensemble averaging (_F[$]) t 
at time t is understood as averaging over the field probabilities pt[$]. This is not 
thermal averaging since we are out of equilibrium. 

The situation we have in mind is one in which, for early times when the available 
space permits many domains, 

(p(x)) t = 0. (5) 

i.e. an equal likelihood of a kink or an antikink occuring in an infinitesimal length, 
compatible with an initially disordered state. However, the density correlation 
functions 

C(x;t) = (p(x)p(0)) t (6) 

will be non-zero, as will 

nit) = (p(x)) t . (7) 

Let us now turn to the complex field </>(x, t) and its more complicated vortices. 
As before we assume that it is possible in principle to calculate the probability 
that 0(x, t) takes the value $(x) at time t f To generalise our observations about 
zeroes on the line to vortices, we follow Halperiol^l in defining the topological line 
density p(r) by 

p(r) = T,Jds^[r-R n (s)). (8) 

In (2.1) ds is the incremental length along the line of zeroes R n (s) (n=l,2,.. .) and 
is a unit vector pointing in the direction which corresponds to positive winding 
number. 

It follows that, in terms of the zeroes of 3>(r), p«(r) can be written as 

Pi (r) = S^ir^dj^d^r), (9) 

where 5 2 [$(r)] = 5[$i(r)]5[$ 2 (r)]. The coefficient of the 5-function in Eq.g is the 
Jacobian of the more complicated transformation from line zeroes to field zeroes. 
We shall also need the total line density p(r), the counterpart of p(x) of Eq||, 

Pi(r) = 5 2 [$(r)] \e ijk d^ x {r)d k ^ 2 (r) | . (10) 

•t Throughout, it will be convenient to decompose $ into real and imaginary parts as $ = -^{^1+1^2) 
(and <f> accordingly). This is because we wish to track the held as it falls from the unstable ground- 
state hump at the centre of the potential to the ground-state manifold in Cartesian field space. 



As before, ensemble averaging (F[Q]) t at time t means averaging over the field 
probabilities p t [$] . Again we assume 

(PiOO>t = 0. (11) 

i. e. an equal likelihood of a string or an antistring passing through an infinitesimal 
street. However, 

n(t) = (pi(r)) t > (12) 

and measures the total string density in the direction i, without regard to string 
orientation. The isotropy of the initial state guarantees that n(t) is independent of 
the direction i. Further, the line density correlation functions 

Qjfrt) = { Pi (v) Pj (0)) t (13) 

will be non-zero, and give information on the persistence length of strings. 

In practice, our ability to construct pt[$] over the whole timescale t > to from 
initial quantum fluctuations to late time classicality is severely limited. It is con- 
venient to divide time into four intervals, in each of which we adopt a different 
approach. If m sets the mass-scale, there is an initial period to <t <t{ — 0(m _1 ), 
before which the field is able to respond to the quench, however rapidly it is imple- 
mented, and which we can largely ignore. Assuming weak coupling, of which more 
later, this is followed by an interval t{ < t < tf in which, provided the quench is 
sufficiently rapid, domains in field phase form, grow, and then cease growing as the 
field magnitude begins to approach the S 1 groundstates. Vortices will appear and 
be driven apart by the coalescence of these domains. In the third period the field 
magnitude relaxes to the ground state values i.e. the vortices freeze in. Finally, in 
the last period, the vortices behave semiclassically. A phenomenological description 
of this period based on time-dependent Landau- Ginzburg theory has been given by 
Liu and Mazenko0, in which they fill in and extend the work of Halperin, and from 
whose calculations we borrow. See also Brayli3. 

Our immediate aim is limited to discussing the second period U < t < tf on 
the (as yet unproven) assumption that the distribution of relevant zeroes at time 
tf is left largely unchanged by their freezing in. This distribution of vortices can 
then be taken as initial data for the final evolution of the network. Although we 
have long term plans for a semi-analytical linkage of this final stage to the initial 
microscopic dynamics, at the very least we can compare the data to that used in 
numerical simulations. 

3. A Gaussian Model for Vortex Distribution 

We have yet to specify the nature of the quench but it is already apparent that, if 
we are to make further progress, additional approximations are necessary. We return 
to our one-dimensional example. The most important concerns the independence, 



or approximate independence, of the fields and their derivatives at the same point. 
If that holds then n(t) of Eq.[7| separates as 



n{t) n (6[*(z)]) t (mx)\) t . (14) 

and we can estimate, or bound, each factor separately. In fact, we can perform 
explicit calculations if p t [$] is Gaussian, for which the approximate equality in Eq.|IJ] 
becomes exact, and this we shall assume henceforth. That this is not a frivolous 
exercise in solving what we can solve but a representation of reasonable dynamics 
will be shown, in part, later, for the 'slow-roll' dynamics that we shall adopt. 

Specifically, suppose that (for the one-dimensional case) $ is a Gaussian field 
for which 

= = ($(x)$'(x)) t , (15) 

and that 

($(x)$(y)) t = W(\x-y\;t). (16) 

All other connected correlation functions are taken to be zero. Then all ensemble 
averages are given in terms of W(r; t) which, from the closed timepath integral 
formalism, can be equally understood as the equal-time Wightman function, 

W{\x-y\;t) = (<l>{x,t)<f>{y,t)) (17) 

with the given initial conditions. In our case, where we shall assume thermal equi- 
librium initially, this is the usual thermal Wightman function. It is straightforward 
to see that 



n(t) = - 

IT 



W"(0;t) 



'11 



W(0;t) 

where primes on W denote differentiation with respect to x. The numerator in Eq.|i~8| 
is the second factor in Eq.[7|. The denominator is the 5-function term, as follows 
directly on writing (5[<&(a;)])t as / da(e ta ^^)t. On using the same exponentiation it 
takes only a little manipulation to cast the correlation function C(x;t) of Eq.||| in 
the form 

c (xi( ) = a*M, (i 9 ) 

where 

Kr,t)= . 7™ . (20) 

The extension to 0(2) line zeroes is messy, but leads to no surprises. Specifically, 
suppose that 

(* a (r)) t = = ($ a (r)^$ b (r)) t , (21) 

and, further, that 

($ a (r)$ 6 (r')) t = W ab (\r - r'|; t) = 5 ab W(\r - r'\;t), (22) 



is diagonal. As before, all other connected correlation functions are taken to be 
zero. 

The density calculation proceeds as before. It follows^" that 
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W"(Q;t) 



W(Q;t) 



(23) 



The line density-line density correlation functions have a much more complicated 
realisationS, albeit still in terms of h(r; t) of Eq.^0], but we shall not consider them 
here in any detail. 

The folly of counting all zeroes is now apparent in the ultraviolet divergence 
of W(r; t) at r = in all dimensions. None of the expressions given above is 
well-defined. To identify which zeroes will turn into our vortex network requires 
coarse-graining, determined by the dynamics. 



4. Gaussian Dynamics and its Coarsening 



It is not difficult to justify our adoption of Gaussian field fluctuations for the 
early period of vortex production. We have already assumed that the initial condi- 
tions correspond to a disordered state. In the absence of any compelling evidence 
to the contrary we achieve this by adopting thermal equilibrium at a temperature 
T higher than the critical temperature T c for t < to- The potential describing this 
state behaves, near the origin in field space, as 

Vffl = \m\T)<Pl (24) 

for an effective mass m(T), with m 2 (T) > . For the sake of argument we take 
V(4>) to be exactly as in Eq.^. The resulting Gaussian field distribution will not 
be seriously disturbed by weak coupling. Further, as we shall see later, initial 
conditions generally give slowly varying behaviour in the correlation function, in 
contrast to the rapid variation due to the subsequent dynamics, and calculations 
are insensitive to them. 

In order to have as simple a change of phase as possible, we assume an idealised 
quench, described by giving a time-dependence to the effective mass m(T) = m(t) 
so that, at t = to, m 2 (t) changes sign everywhere. Given that the resulting theory 
displays symmetry-breaking, this change in sign in m 2 (t) can be interpreted as 
due to a reduction in temperature. Even more, we assume that, for t > to, m 2 (t) 
takes the negative value m 2 (t) = —M 2 < immediately, where — M 2 is the mass 
parameter of the (cold) relativistic Lagrangian. That is, the potential at the origin 
has been instantaneously inverted, breaking the global 0(2) symmetry, interpreted 
as a temperature quench from high-T to, effectively, zero temperature. If, as we 
shall further assume, the A|0| 4 field coupling is very weak then, for times M(t — to) < 
ln(l/A), the 0-field, falling down the hill away from the metastable vacuum, will 
not yet have experienced the upturn of the potential, before the point of inflection. 



Thus, for these small time intervals, pt[3>] is Gaussian, as required. Henceforth, we 
take t = 0. 

For such a weakly coupled field the onset of the phase transition at time t — 
is characterised by the instabilities of long wavelength fluctuations permitting the 
growth of correlations. Although the initial value of (</>) over any volume is zero, we 
anticipate that the resulting evolution will lead to domains of constant (0) phase, 
whose boundaries will trap vortices. 

Consider small amplitude fluctuations of <p a , at the top of the parabolic potential 
hill. Long wavelength fluctuations, for which |k| 2 < M 2 , begin to grow exponen- 
tially. If their growth rate Q(k) = y 1 M 2 — |k| 2 is much slower than the rate of change 
of the environment which is causing the quench, then those long wavelength modes 
are unable to track the quench. Unsurprisingly, the time-scale at which domains 
appear in this instantaneous quench is £j = 0(M _1 ). As long as the time taken to 
implement the quench is comparable to and much less than tj = 0(M~ 1 ln(l/X)) 
the approximation is relevant. 

This picture essentially tells us how to coarsegrain the ultraviolet zeroes to iden- 
tify the vortices that will form the network. Firstly, we note that if $ is Gaussian, 
then so is the coarsegrained field on scale L, 

$ L (r) = J c/V/(|r-r'|)$(r'), (25) 

where I(r) is an indicator (window) function, normalised to unity, which falls off 
rapidly for r > L. The only change is that Eq.^2| is now replaced by 

($L,a(r)$L, b (r')> t = W L , ab {\r - r'\;t) = 5 ab W L (\r - r'\;t), (26) 

where Wl = J J IWI is now cut off at distance scale L. Wl(0] t), its derivatives, and 
all relevant quantities constructed from Wl are ultraviolet finite. The distribution 
of zeroes, or line zeroes, of $^ is given in terms of Wl as in the previous section. 

Choosing L = M~ l solves all our problems simultaneously. At wavelengths 
k~ l < L i.e. k > M) the field fluctuations are oscillatory, with time scales 0(M _1 ). 
Only those long wavelength fluctuations with k^ 1 > L have the steady exponential 
growth that can lead to the field migrating on larger scales to its groundstates. Fur- 
ther, as the field settles to its groundstates, the typical vortex thickness is 0(M _1 ), 
and we only wish to attribute one zero to each vortex. By taking L = M _1 we are 
choosing not to count zeroes within a string, apart from the central core. 

We are now in a position to evaluate pt[$], or rather Wl{t; t), which we now write 
as WM{f",t) for t > 0, and calculate the defect density accordingly. This situation 
of inverted harmonic oscillators was studied many years ago by Guth and PiUl and 
Weinberg and WuS. In the context of domain formation, we refer to the recent 
work of Boyanovsky et al.i, and our own!. For weak coupling we recover what 
would have been our first naive guess for the coarse-grained correlation function 
(<t>L,a{ r ,t)(f)L,b{ r ',t)) based on the growth of the unstable modes a (k, t) ~ e n ( fc )* 
alone, 

W M (r;t)~ [ fke ikr e 2n{k)t , (27) 

J|k|<M 



provided we have not quenched from too close to the transition §and t > M _1 
(preferably t 3> M _1 ). The reason why Wm{j"i^) ls approximately independent of 
the initial conditions is that the integral at time t is dominated by a peak in the 
integrand & 2 e 2n ( fc )* at k around k c , where tk\ « M. Any temperature dependence 
only occurs in slowly varying factors (like coth(3u)/2) and can be ignored. 



5. Densities and Correlations 



In the first instance we understand the dominance of wavevectors at k 2 c = M/t 
in the integrand as defining a length scale 



£(t)=0{y/i/M), (2i 



once Mt > 1, over which the independently varying fields (f) a are correlated in 
magnitude. However, in our introductory comments we anticipated that string 
distributions would be determined by the length over which the phase of the field is 
correlated. For this simple model, with only a single length, £(£) has to serve both 
purposes. 



To see how this happens, we take = 2Jt/M. To calculate the number 
density of vortices at early times we insert the expression Eq.|27] for Wm into the 
equations derived earlier, to find 

7r£(t) 2 

permitting us, from our introductory comments, to interpret £(t) as a correlation 
length for phases. We note that the dependence on time t of both the density and 
density correlations is only through the correlation length We have a scaling 
solution in which, as the domains of coherent field form and expand, the interstring 
distance grows accordingly. Since the only way the defect density can decrease 
without the background space-time expanding is by defect-antidefect annihilation, 
we deduce that the coalescence of domains proceeds by the annihilation of small 
loops of string, preserving roughly one string zero per coherence long held 

belief for whatever mechanism. 

As well as determining the string density, we are also looking for signs of anti- 
correlation, which enables us to determine the persistence length of strings in the 
network (i.e. how bendy they are). The bendier, the more string that will occur in 
small loops. This is important, since the conventional string model for large scale 
structure formation in the universe requires a certain amount of infinite string. In 
fact, there may be no need to calculate the line density correlation functions CV, to 
appreciate that there is a higher fraction of string in loops than would have beeji an- 
ticipated from the early simulations which laid out a regular domain structured. On 



§In which case initial fluctuations are so large that the field magnitude has already sampled the 
ground states and a slow roll is inappropriate. 



inspection, the broad peaking in wavelength I = k~ x about k = k c in the integrand 
of Eq.^ is understood as indicating that the domains, both in field magnitude and 
phase, of characteristic linear dimension have a substantial variation A£(i) in 
size. Simple calculation shows that 

m ~ k c ~2- (30) 

independent of time. 

To make use of Eq.|30] requires two further assumptions. The first, which can 
be checked using the methods of Ref.7, is that the domains growing as stop 
expanding almost immediately as soon as the field magnitudes probe the spinodal 
region and the instability switches off. The second is that, during this time and 
the time it takes for the field magnitude to freeze on the S 1 groundstates, A£/£ 
of Eq.^0] describes (albeit approximately) the variation in domain size over which 
field phase is correlated. In numerical simulation of string networks the inclusion of 
variance in the size of field phase domains showdll that, the greater the variance, 
the more string is in small loops. The variance of Eq.|30| suggests much more string 
(e.g. twice or more) in small loops. However, it is not easy to marry the somewhat 
different distributions of this simulation to ours. We can say no more than there 
seems to be some infinite string and, in an early universe context, some is enough. 

To do a little better we need to evaluate density correlation functions. For Wm 
above we find that, in units of n 2 (t), the anticorrelation is large. In particular, a 
calculation of C33 for separation r in the 1-2 plane gives 

C33M) = ~h 2 {r;t) 



r dr 

At) 



-1 + 



e(t) 



(31) 



when r < £(£). Although it is difficult to be precise, this again suggests a significant 
amount of string in small loops. For r ^> £, C33 falls off exponentially as 

C 33 (r;t) = 0(e- 2r2 /^). (32) 

showing that £ indeed sets the scale at which strings see one another. Unfortunately, 
it has not yet proved possible to turn expressions for the Cy into statements about 
self-avoidance, fractal dimension, or whatever is required to understand the length 
distribution of the resulting string network. 

However, we can use the correlation functions to determine the variance in vortex 
winding number Ng through an open surface S in the 1-2 plane, which we take to 
be a disc of radius R. There is a complication in that, as can be seen from Eq.|8], p 3 
counts zeroes weighted by the cosine of the angle with which they pierce S, whereas 
winding number counts zeroes in S without weighting. The problem is therefore 



essentially a two-dimensional problem. Using topological charge conservation across 
the whole plane, it is not difficult to see that 



(AN) 2 S = 0(n 2 e(t)R) = 0(R/at))- 



(33) 



All that is required is short range zero-density correlation functions. Since (AN)s 
is (Aa)s/27r, where (Aa)s is the variance in the field phase around the perimeter 
dS of S, this means in turn that 



This is naturally interpreted as a random walk in phase along the perimeter with 
average step length within which the phase is correlated. Calculations will be 
given elsewherea. 

Since only the boundary region of S gives a contribution to the variance, we 
would get the same result if we were to quench in a annulus. It has already been 
proposed! that quenches in an annulus be performed for superfluid A He, for which 
(Aa)s can be identified with a supercurrent. 
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